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High-energy collisions of two nucleons on two nucleons are studied in the BFKL approach in the 
leading approximation in a s N c . Diagrams with redistribution of colour are considered. It is found 
that intermediate BKP states consisting of 4 reggeized gluons give a contribution which may be leading 
in deuteron-deuteron scattering and thus experimentally observable. 
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Collisions of two heavy nuclei have long occupied a prominent part of the experimental and theoretical 
studies in strong interactions at high energies. Unlike the case of DIS theoretical analysis of these 
\ processes turned out to be quite complicated. The most advanced calculations have been made in the 
framework of the Colour Glass Condensate approach, where they heavily rely on numerical Monte- 
Carlo method on the lattice to study evolution of the classical gluonic field in the course of collisions 
pJ,GJ[3]. In comparison, analytic methods applied to collision of heavy nuclei have only given modest 
and approximate results [H [5j |6] . So to understand the problem it seems natural to start not from this 
general case, but from the simplest generalization of the well-know results for DIS to the collision of 
two nucleons with two nucleons. The immediate physical application is of course the case of deuteron- 
deuteron collisions. An alternative view is to consider this as a particular contribution to the amplitude 
for the collision of two heavy nuclei A and B, which generally contains contributions from collisions 
of any number of nucleons in the projectile nucleus on any number of nucleons in the target. For 
large atomic numbers A and B the scattering amplitude is effectively unitarized by eikonalization of 
its connected part E(b) of the corresponding diagrams: 

AAB = 2isJd 2 b(l-e~ E ^), (1) 

where b is the impact parameter of the collision and s is the cm energy squared for a pair of colliding 
nucleons. Our study then refers to the part of the eikonal function coming from the collision of two 
nucleons from the projectile nucleus on two nucleons from the target nucleus. Due to eikonalization 
the total nucleus-nucleus cross-section for heavy nuclei does not practically depend on a concrete value 
of the eikonal function. The latter is then large and its changes even by several times do not mean 
much for the total cross-section, which remains essentially geometrical. However even for heavy nuclei 
it may have influence on some specific processes, like exclusive rare production, where it may change 
the absorptive factor, and of course on the inclusive cros-sections. 

It is instructive to see relative orders of magnitude for various contributions to the eikonal function 
E(b). Naturally the leading order in both coupling constant a s and number of colours N c is given 
by the double gluon exchange, Fig. [Lj a. For simplicity, instead of nucleons we shall consider quarks 
as elementary scattering centers in the nuclei, invoking the colour neutrality by projecting onto the 
colourless t-channels. 
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Figure 1: Lowest order diagrams for the eikonal: two gluon exchange (a) and with transitions to 2 (b) 
or 3 (c,d) intermediate gluons 

Note that the total diagram for 2 by 2 scattering involves also the product of two such diagrams. 
However the eikonal function will include only one of them. If we forget about couplings inside the 
scattering centers, its contribution does not depend on a s and is proportional to N 2 . Other connected 
diagrams may be classified by the minimal number of exchanged gluons n m i n . We shall discuss their 
order of magnitude related to the double gluon exchange, Fig. Q] a. The diagram in Fig. [T] b with 
n-min = 2 has this relative order a 2 N 2 . The diagram in Fig(Hc with n m i n = 2 has the relative order 
agN^. Diagrams d and e with n m i n = 3 have order a^N 3 . Finally typical diagrams with n m i n = 4 and 
redistribution of colour, which makes them connected, are shown in Fig. [2] and have the same order 
as Fig. [TJi, so that their relative order is unity. To finally estimate the weights of these diagrams one 
has to take into account that diagrams with n m i n = 2 and 3 involve one or two intermediate rapidities 
at which the initial and final 4 gluons fuse into the intermediate ones and so are proportional to y or 
y 2 where y is the overall rapidity. So their final order will be a 2 N 2 y, a^N^y 2 and a^N^y 2 . In the 
BFKL kinematics one assumes a s N c y ~ 1, a s N c « 1 so that orders of diagrams Fig. QJ, c, d and 
e become a s N c , a 2 N 2 and a s N c respectively. This shows that apart from the double gluon exchange 
the dominant contribution comes from the diagrams with colour redisitribution, Fig. [2j Of course 
this result has been obtained in the lowest order, but it remains valid also in higher orders when the 
leading order will be multiplied by powers of a s N c y ~ 1. 

In this study we shall consider the connected part of the forward scattering ampitude of two 
projectile centers (quarks) on two target centers (also quarks). As we have demonstrated, in the 
lowest order it is just the simple rearrangement amplitude shown in Fig. [2] in two different forms: 
one (a) symmetric in projectiles and targets and another (b) showing the intermediate states in the 
s-channel. In fact both in Fig. [2] a and b one should also take into account all diagrams with crossed 
vertical lines (16 diagrams in all). We would like to study all possible inclusions of interactions between 
the gluons. They are realized by the BFKL interactions Vij between gluond i and j [7J. We shall work 
in the leading approximation in the number of colours N c . As compared to the disconnected diagram 
with two pomeron exchanges of the leading order in N c , all contributions studied in the following will 
be of the order 1/N 2 and so subdominant in the large N c limit. However in the eikonal function they 
will contain an extra nuclear factor of the relative order A 2 / 3 for collisions of two nuclei of the same 
atomic number A and so in fact may be of the same order or even greater than the leading contribution 
in N c . 

We can separate all contributions into three classes. First, interactions may occur only between 
gluons attached to the same projectile or target and so being in the vacuum state. Such interactions 
will convert pairs of gluons attached to the projectiles and targets in Fig. [2] together with their 
crossing into fully developed pomeron Green functions, so that the quarks representing our projectiles 
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Figure 2: Lowest order diagram with a redisiribution of color (diagrams with crossing gluon lines 
should be added) 
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Figure 3: Diagrams for the direct transition of pomerons with the redistribution of colour 

and targets will convert to pomerons, see Fig. [3j For heavy nuclei reggeized gluon splitting will further 
convert simple pomerons into fan diagrams made of pomerons, which are summed into the solution of 
the Balitski-Kovchegov (BK) equation, BK wave functions [3 [9]. 

Second class of the diagrams are formed form those of the first class with one interaction between 
the gluons V23 or V14, which does not connect any pair of gluons in the vacuum state (Fig. 

Third class of the diagrams are those where from each side there appear interactions which connect 
vacuum pairs of gluons with the so-called BKP states |10t 111] between them. (Fig. [5|). 

Before actually calculating all these contributions, in Section 2 we develop a multirapidity formal- 
ism which makes it easier to sew four pomerond into a single amplitude. 

Let the momenta of the nucleons in the projectile be k and those in the target be I. In the cm. 
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Figure 4: Diagrams with the redistribution of colour and one interaction beiween the pomerons of the 
projectile and target 
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Figure 5: Diagrams with the intermediate BKP states 

system of a pair of nucleons from the projectile and target we have 

k+ = I—, k- = k± = l + = l± = 0, s = 2k + l- 

The amplitude for nucleus-nucleus scattering can be separated into its high-energy part H and two 
nuclear factors (see Appendix 1.). The imaginary part of H in its turn can be presented in the form 

ImH = -(27r) 2 5(n + )5(q-)4s 2 N?D, (2) 

where k and q are the momenta transferred ito the projectile and target nuclei respectively with 
k_ = k± = q + = q± = and D is a real function which is the sum of connected diagrams for the 
forward scattering of two-nucleons on two nucleons written as integrals in rapidity and transverse 
momentum space. At fixed impact parameter b the cross-section for the part of nucleus-nucleus 
scattering coming from interactions of two nucleons from the projectile on two nucleons from target 
is related to D as 

af B (b) = -2N 2 C D J dH A Tl(b A )Tl(b - b A ) (3) 
and for deuteron-deuteron scattering 

(see Appendix 1.). 

Note in conclusion that the contribution of simplest diagrams with colour rearrangements have 
been considered in literature in relation to correlations between a pair of produced gluons |12| . We 
postpone discussion of the inclusive gluon production to future publications, since in our BFKL- 
Bartels approach it involves the study of possible cuts of the forward scattering amplitudes, which 
for 2 by 2 scattering requires special attention to the well-known AGK cancellations. Meanwhile we 
have to stress that already on the level of the total cross-section, apart from the pomeron (or BK) 
wave function, our amplitudes involve more complicated objects made of 4 reggeized gluons (the BKP 
states). In the Colour Glass Condensate approach similar complicated structures begin to appear only 
on the level of double inclusive cross-sections. 



2 Pomeron in the multirapidity formalism 

In this section we introduce a formalism in which each pomeron has its own rapidity, which corresponds 
to the standard Feynman diagram technique. In this formalism construction of amplitudes with 
colour rearrangment becomes much simpler. Recall that in the BFKL approach the amplitude ("wave 
function") P(y) at a given rapidity y is obtained from its value at y = by the transformation 

P(y) = e-y H P(0). (5) 



Here the BFKL Hamiltonian is H = H^ ' + V. The unperturbed Hamiltonian is a sum of Regge 
trajectories u(k) with a minus sign. For gluons 1 and 2 H-j^ = ~ ^2 or, in the momentum 
representation, 

< k' x ,k' 2 \H^\k u k 2 >= -(2ir) 4 5 2 (k 1 - k[)5 2 (k 2 - j^)(w(ifei) + u(k 2 )). (6) 

The potential energy V is given by the pair BFKL interaction. Symmetrizing in initial and final 
reggeons we have 

V12 = {T 1 T 2 ){2^) 2 5 2 {k 1 + k 2 - k[ - tel**, fc), (7) 

where 

.2 , l2 U 2 _i_ 1,2 U 2 



J? / fcl ^2 + ^2 ^1 

2nk 1 k 2 k' 1 k' 2 \ (h - k[) 



k) = ^ ^ 121 - (A* + fc 2 ) 2 (8) 



and Ti and T 2 are colours of the two gluons. 

Evolution law ([5]) mimics the standard evolution in time provided one makes the substitution 

t -> y -»• y > 0. (9) 

By definition P(y) = at y < 0. We shall obtain our diagrammatic technique studying evolution in 
time. Evolution in rapidity will be obtained making the analytic continuation ([9]) at y > in final 
formulas. 

Passing to the multirapidity formalism we shall additionally characterize each reggeon by its "en- 
ergy" e. We shall take the reggeon propagator as a function of e in the form 

A(e,k) = t~t . (10) 

Then as a function of time 

A (t , k) = if^—±^ = me ^K (11) 

Taking it — > y at t > and y > we get the desired reggeon propagator 

A(y,k) = e( y y»W. (12) 

The interaction between reggeons will be given by ([7]) with factor (— i) and with an additional 
factor responsible for conservation of the total energy 

V X2 = -i{T 1 T 2 ){2n) 3 S(e 1 + e 2 - e[ - + k 2 -k' 1 - k^v^^lM)- (13) 

These Feynman rules allow to construct amplitudes for the interaction of any number of reggeons 
provided this number does not change. 

Our first task is to check that for a pair of reggeons the sum of all thus consructed Feynman 
diagrams is equivalent to the standard BFKL equation. 

The Green function for a pair of interacting reggeons Q{^\k' x ^ 2 k 2 \e\k\^e 2 k 2 \ which is illustrated 
in Fig. [H obeys an equation which is graphically shown in Fig. 

Separating from Q the 5-functions which corresponds to energy-momentum conservation 

^(e / 1 /c / 1 ,e / 2 A: 2 |eiA;i,e 2 A; 2 ) = (2n) 3 5(e' 1 + e' 2 - ei - e 2 )5 2 (k[ + k' 2 -k\ - k 2 )GEK(exk[\exki), (14) 

where E = t\ + e 2 = e' x + e 2 and K = k% + k 2 = k[ + k' 2 . we find a Bethe-Salpeter type equation for 
Qek 

GEK^k'^h) = (27r) 3 8(e' 1 -e 1 )S 2 (k' 1 -k 1 )A 1 A 2 + A' 1 A' 2 J ^^V(k{XWlX)SEK^X\^i), 

(15) 




Figure 6: The pomeron Green function Q 





Figure 7: Equation for the pomeron Green function Q 



where we denote Ai = i/{e\ +u>i), A[ = i/{e'i and so on and all w's are assumed to have a small 

positive imaginary part. 

In this equation the interaction is energy independent, as in theories with non-relativistic potentials. 
So it can be easily transformed into a Schroedinger-like equation. Indeed if we present 



QEKie'A^h) = (27r) J <J(ei - ei )5 2 0i - k 1 )A 1 A 2 + A' 1 A' 2 T EK A 1 A 2 



(16) 



then Tek (e'ik[\eiki) does not depend on initial nor final energies e±, e 2 , e[ and e' 2 . The equation for 
T takes the form 



f EK (k' 1 \k 1 ) = V(k' 1 ,k , 2 \k 1 ,k 2 ) + i 



d 2 k'{ V(k[,k 2 \k'{, k%)T EK (k![\kx 



(17) 



Recalling that V 
operatorial form 

where 



' (2tt) 2 E + u'{ + J 2 ' + iO 
-iV, where V = (TiT 2 )v, and putting T = —iT we rewrite Eq. (|17p in the 

T EK = V + VR EK T EK (18) 



Rek = {E- # (0) + iO)" 1 



(19) 



is the resolvent of the unperturbed Hamiltonian. Eq. (|18j) is the standard equation for the T-matrix. 
The standard BFKL Green function is defined by T E k as 

G EK = Rek + RekT ek Rek = (E- H)' 1 (20) 

The Bethe-Salpeter Green function Q is expressed via the Schroedinger one G, as 

QEKiA^iki) = {2^f5{e[-e 1 )5 2 {k[-k 1 )A 1 A 2 + i{2^f5 2 {k[-k 1 )[E + u l +oj 2 ]A' 1 A' 2 A 1 A 2 

- iA[A' 2 [E + u[ + oj'^GEKik'^k^lE + ^ + u^A^. (21) 

Now we pass to the BFKL function proper. Integrating the Green function Q E K( e ik'x\eiki) with the 
impact factor p E K(k' 1 )/(k' 1 k' 2 ) we obtain a function which describes the pomeron in the multirapidity 
formalism 

r df' ii 2 h' 

VEK^h) = J £ 3 ^ PEK(K)SEK&K\*lkl). (22) 
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Figure 8: Equation for the pomeron wave function V 



The equation for it easily follows from Fig. [8] 



V E K{exk 1 )=vf K {e 1 k x )-i^ 2 J d ^^V{k 1 M\k[A)'PEK{e' 1 k' 1 ) 



(2vr) 3 



(23) 



where e 2 = E — ei, k 2 = K — k\ and oj 2 = u(k 2 ). If we introduce the "amputated" wave function by 



V E K(eih) = F EK {e u fei)AiA 2 
then £ does not depend on energy ei and satifies 

d 2 k[ vfaMK^FEK&K) 



F EK (h) = FjVfa) + J 



So if we define a new wave function 



P EK (ki) 



(2vr) 2 E + oj[ + oj' 2 

F EK (ki) 



E + OJi + U2 



it will obey the equation 

(E + uj 1 +u 2 )P EK (k 1 ) = F^ K (k 1 )+ J ^ I V(k l Mk'i,k' 2 )P EK {k' l 



(24) 



(25) 



(26) 



(27) 



which is the standard BFKL equation (for the pomeron or gluon depending on the value of (TiT 2 ) in 
V). 

Note that BS function V turns out to be related to P as 



V EK (t\k\ 



E + + u 2 
(ei + wi)(e 2 + oj 2 ) 



P EK (ki 



(28) 



It depends on the individual energies. It is remarkable that although for the pomeron P is infrared 
safe, the coresponding BS function V does not look so. 

3 Four pomerons sewed with the redistribution of colour 

Now we pass to the simplest amplitude D\ in which two pomerons from the projectile (reggeon pairs 
(1,2) and (3,4)) are directly coupled to two pomerons in the target (reggeon pairs (13) and (24)) 
illustrated in Fig. [3l All pomerons will have their total momenta equal to zero, which will be tacitly 
assumed in the following. The pomerons from the projectile have their total energies E\ 2 and £^4 
and those from the target energies E\% and £24- Energy conservation requires £"12 +£34 = £13 + £24- 
Correspondingly the amplitude D\{E\ 2 , £34, £13, £24) will include factor 2tt5(Ei 2 + £34 — £13 — £24) 
which will be suppressed in the following. As seen from Fig. [3l the internal integrations will include 
a single transverse momentum q and a single pomeron energy, say, e\. The remaining energies will be 
expressed via ei as follows 



e 2 — £12 — ei, £3 — £43 — ei, £4 — £24 — e 2 — £24 — £12 + ei- 



(29) 



In terms of the amputated pomeron function F the amplitude will be given by 

n (17 w T? F ^ f deid 2 q F El2 (q)F E . M (q)F El3 (q)F E24 (q) 
D1 (E 12 ,E U ,E 13 ,E U ) - J (27r)3 (ei + w (,))( e2 +a;( g ))(e3+a;(9))(e 4 + a,(,))' (30) 

where all w's are assumed to have a small positive imaginary part. 

Since F's do not depend on e\ we can do integration on ei explicitly. We have an integral 

da 1 1 l l 



2vr ei + + iO £? i2 - e x + + iO £7 13 - ei + + iO £ 2 4 - #12 + ei + w(g) + iO ' 
i r 1 1 1 



£ 13 - £12 I E 12 + 2oj{q) + iO E 24 + 2w(q) + iO £13 + 2w(g) + iO £ 3 4 + 2w(g) + iO }" ^ 

Now we express our amputated pomerons i* 1 via full ones P using Eq. ([260 . In terms of P the 
amplitude will be given by 

D l {E 12 ,E^,E l ^E 2A ) = -iJ —^^E 12 + E 3 4 + M<l))PE 12 (q)PE 3 M p Eis(Q)PE 2 M- (32) 

Immediately the question of its infrared safeness arises, since the integrand contains the reggeon 
trajectory oj{q). 

Next we study evolution in rapidity. 

DM = J dEudE^dE^ ^^ + ^ _ E ^ _ E24)e -v(E 12+ E. iiDi{E ^ E ^ ^ ^ = 
■ f d 2 q (a ( ^ <9x f dE 12 dE u dE 13 dE 2U 

~ l i J2^( 4u{q) -an) J (2^ 27T6{El2 + Eu ~ El3 ~ Eu) 

^y^+E^p^ {q) p Eii {q) p E2i {g) (33) 

We first consider evolution of the integrand of the momentum integral in time: 
m = J dE 12 dEudE l3 dE 2 ^ §{Ei2 + ^ _ Ei3 _ £24)e -, t (% + %) P£i3( ^ £34(g)P£i3(?)?%(g) 

= J dE 12 dEudE 13 dE 2i dTeiT(gia+gM - gl ,-^) e -it( gl3+ ^) P£ , ia {q) p Em {q) p Ei3 {q) p E24 (g) . (34) 

J dEP E {q)e- itE = 6(t)P(t,q) (35) 



We use 



to obtain 



Z(t) = J dr9(t - r)P 2 (t - r, q)9(r)P 2 (r, q) = 9{t) J* drP 2 (t - r, q)P 2 (r, q). (36) 
Analytic continuation to rapidity gives 

Z(y) = -i I dy'9(y - y')P 2 {y - y', q)9(y')P 2 (y' , q) = -i9(y) f dy'P 2 {y - y\ q)P 2 (y', q), (37) 



so that finally 

Dl (y) = -9(y) J ^(4u(q)-l) jT dy> P 2 (y - y> \q)P 2 {y' , q). (38) 

To check the correctness of the transition from time to rapidity we study a simple model, in which 
P E is given by a BFKL pole 

P E = , (39) 

E + a + iO V ' 



so that P(y) = 6(y)e ay . In this case it is easy to find evolution of Z(y) in rapidity by explicit analytic 
continuation. Let for complex z 

z{z) = J dEdEfodEu e - zE p Ei2 {q) p E _ Ei2 {q) p Ei3 ( q )P E _ Ei3 (g). (40) 



We have 



dE> rdE> 



So we find 



2vr J 2vr (E> + a + iO)(E - E> + a + iO) E + 2a + i0 

But this integral only exists for pure imaginary z = it when it is given 

Z (t) = - [ —e~ itE - ^ = +ie(t)4-e- itE \ E =-2a = te 2aU -> -ize 2az . (42) 

w 7 2vr (£ + 2a-M0) 2 w aE 1 v ; 

On the other hand, from (13711 we find 



Z(y) = -i0(y) f V dy'e 2a{ - y - yl) e 2ay ' = -iye 2ay (43) 
Jo 

in accordance with ([U]). This confirms factor (— i) in ([37|) and so ([38]) , 

4 Single interaction between the projectile and target pomerons 

Now we consider diagrams with a single interaction between the projectile and target which cannot 
be included into the pomerons, that is V23 or Vn, see Figs. HJi and b. 

We start with the diagram shown in Fig. SJi. As before we suppress the energy conservation factor 
2ir5(Ei2 + £34 — £13 — £^24)- The diagram contains two loops and so integrations over ei, €4,(71 and 
q 4 . Energy-momenta of the gluons 1,2,3 and 4 before the interaction are 

(ei,qi), (E 12 - e 1 ,-q 1 ), (E 34 - e 4 , -q 4 ), (64,54). (44) 

After the interaction gluon 2 and 3 have their energy-momenta (E24 — 64) and (-E31 — e\, —qi) respec- 
tively. In terms of amputated pomerons F the contribution from Fig. [4] a is 



E>2a = i 



de 1 de 4 d 2 qid 2 q 4 F El2 (q 1 )F E34 (q 4 )F El3 (q 1 )F E24 (q 4 )V(-qi, -q 4 \ - q 4 , -qi) 



(27r) 6 (ei + ui)(E 12 - e\ + ui)(E 13 - ei + wi)(e 4 + W4X-E34 - e 4 + w 4 )(£ 2 4 - £4 + W4) ' 

(45) 

The prefactor includes — i from the definition of V" and (—1) from 6 propagators. 
Integration over energies factorizes into two integrals: 



de± 



2vr(ei + wi + i0)(^i 2 - ei + wj + iO)(E 13 -e 1 + u 1 + iO) (E 12 + 2wi)(£?i 3 + 2wi 
and a similar integral over €4 which gives 

1 



(46) 



(47) 



(£34 + 2lu 4 )(E 24 + 2w 4 )' 
So we get 

^ _ _ { f d 2 qid 2 q4F El2 (qi)F E . M (q 4 )F El3 (q 1 )F E24 (q 4 )V(-q 1 , -q 4 \ - q 4 , -qi) 
2a J (2tt) 4 (E 12 + 2u; 1 )(E 13 + 2uj 1 )(E u + 2uj 4 )(E2 4 + 2oj 4 ) 1 ' 



Recalling relation ([26]) between the amputated wave function F and pomeron we rewrite this as 
D 2a = -i J 7^34 4 PE 12 (qi)PE 34 (q4)PE 13 (qi)PE 24 (q4)V(-qi,-q4\ - g4,-gi). (49) 

Now we pass to the amplitude corresponding to the diagram in Fig. [4] b. Before the interaction 
gluons 1,2,3 and 4 have the same energy-momenta as before (Eq (|44|) ). However now the momenta 
of gluons 1 and 4 change after the interaction and become (Ei$ — £34 + £4,94) and (E24 — E\2 + ei) 
respectively. 

In terms of F the amplitude is 

D2b = l Wf x 

f de 1 de 4 d 2 q 1 d 2 q 4 F El2 (qi)F E34 (q 4 )F El3 (q 4 )F E24 (q 1 )V(q 1 ,q 4 \q 4 ,q 1 ) 

J (ei + u\){E\2 ~ ei + wi)(-E , 24 - ^12 + ei + wi)(e 4 + w 4 )(-E 34 - e 4 + ^ 4 )(Ei 3 - £34 + e 4 + w 4 ) ' 

(50) 

Again the integrals over energies factorize into two ones: 

j = f dei = 1 

3 J 2vr(ei + wi + iO)(-Ei 2 - e x + wi + i0)(£ 24 - E l2 + ei + wi + iO) (Ei 2 + 2wi)(£ 24 + 2wi) 

(51) 

and a similar integral over £4 

h = i-n=, Itt=, c- (52) 

(£34 + 2^4) (£13 + 2^4) V ' 

We get 

D __. f d 2 qid 2 q 4 F El2 (q 1 )F E34 (q 4 )F El3 (q 4 )F E24 (q 1 )V(q 1 ,q 4 \q 4 ,q 1 ) 
2b 1 J {2^Y{E 12 + 2uj 1 ){E^ + 2uj 4 ){E Z4 + 2uj 4 ){E2 4 + 2uj 1 ) 1 ' 

or in terms of pomerons P 

4 p Ei2 (qi) p E M {qa)Pe 13 {qA)PE 2 A (gi)y(gx, g4[g4, gl)- (54) 

Let us separate the infrared stable and divergent parts in D 2 - For a pair of gluons 1 and 2 the 
BFKL Hamiltonian is 

H = -U)i -uj 2 + (TiT 2 )v 12 , 

where v\ 2 is given by ([8]). In the vacuum t-channel (TiT 2 ) = —N c , so that the infrared stable pomeron 
Hamiltonial is 

Hp = -u>i - oj 2 - N c vi 2 . (55) 

In D2 interaction connects different colour configurations with gluons from the projectile forming 
colourless pairs (1,2) and (3,4) and from the target forming colorless pairs (1,3) and (2,4). The 
transition matrix element of (T2T3) entering the interaction is +N C (see Appendix 2.). So in terms of 
v amplitude I? 2a is given by 

— T^i - p E 12 (qi)PE34 {<&)Pe 13 (gi)PB 24 (g4)Ar c v(-gi, -g 4 | - g 4 , -qi). (56) 

We present 

N c v(-qi, -g 4 | - g 4 , -gi) = - < -gi, -g 4 |-ffp| - g4 - gl > -2wi(27r) 2 (5 2 (gi - g 4 ), (57) 

so that 

f Q\ cfi Q^. 

D 2a = i J (2^)4 4 Pei2 ^ Pe 34 (<te) p E 13 (qi)PE 24 (94) < -qi,-qi\H P \ - g 4 ,-gi > 



1 n 



P d? 

+ 2i J J^^) P ErMPE z MP El MPE 2 M (58) 



(27T) 
and similarly 

/d*^ Q\ ^4 
(2tt) 4 PEi2 ^ PE34 ( g4 ) PE i3 (Q4)Pe 2 4 (gl) < 9l,94|#p|94,9l > 

+ 2i/ ^2^(9)^ 12 (9)^34 (9)^3(9)^4(9)- (59) 

As we see the additional terms in (|58p and (|59p containing w(g) cancel with a similar term in (|32p . so 
that the remaining sum of Z?i and D 2 turns out to be infrared safe. 

So after cancellation of the gluon Regge trajectories the infrared safe contributions from diagrams 
in Figs. [3] and H] are 

Dl = 6{v) i) \l dV 'I ^ P2 ^-^) P V,9) (60) 

and 

D 2 = 26(y) J'dy'J ^0jf < q,q'\H\q' ,q > P(y - y> ,q)P(y - y> ,q')P(y' ,q)P(y' ,q>). (61) 

5 Two interactions between the projectile and target pomerons 

With two interactions between the target and projectile pomerons we can use our old results in [13] 
where coupling of two pomerons to the BKP 4-gluon case was studied (see also Appendix 2.) 

In this case there are transitions both with the redistribution of colour, that is | (12) (34) >— > 
I (13)(24) >, and without this redistribution, that is |(12)(34) >^ | (12)(34) >. 

To continue our line of studies we start with the redistribution of colour. As follows from our 
studies in Appendix 2 in this case between the pomerons from the projectile and target there can 
appear two 4-gluon BKP states 1 1243 > and 1 1324 >. Inserting the 4-gluon BKP Green function G 
between them we find that the projectile and target pomerons will be connected by 



M ( E a) 



-N 2 C (v 13 + v 24 - v 23 - v M ) [Gg 243) + 4 1342) ] (« 12 + v u - v 23 - v u ) , (62) 



where, say, G^ 243 ^ is an operator acting in the 4-gluon space 

1^,(1243), / / / / . 

< 91, 92,93, qA G e l9i, 92,93, 94 > 

satisfying the equation 

(£_tf(1243) )G 1243 = 1 ( g 3) 

with 

4 1 

#(1243) = _ Ui - -g 2 N c (v 12 + v 24 + ^43 + «3l)- (64) 

i=l 1 

Note that can also be presented in terms of the infrared safe BFKL Hamiltonian for the pomeron 
Hp 

M E ] = J ( H PA3 + H P,24 ~ Hp,23 ~ H PM ) [G£ 243) + Gf 2A) ] (h p>12 + H PtU - H P , 23 - H PM ) (65) 



which demonstrates that M E is infrared safe 



(a) . 
E 

The explicit form for the kernel of M is 

4 , , J .;-'/./ 



< 91, 92, 93, 94|i< \q[,q' 2 , q>, q' A >= 2n5* ( £ 9^E *) \ N c / II ^ ^f] 2 ^ ( E *f"E **) 



-1 -I 



< Qi,Q2,Q3,Q4\vi3 + ^24 - v 2 3 - vu\h, k 2 ,k 3 ,k 4 >< h,k 2 , k 3 , A; 4 |G^ 243) + G^ 342) | k[ , k' 2 ,k' 3 ,k' 4 > 

< k[, k' 2 ,k' 3 ,k' 4 \v 12 + t>34 - v 23 - vu\q[,q' 2 ,q 3 ,q 4 >> ( 66 ) 

where for instance 

< QiiQ2,q3,q<l\vi3\ki, k 2 , k 3 ,k 4 >= (27r) 4 (5 2 (<72 - k 2 )5 2 (q i - k 4 )v(qx,q 3 \k\,k 3 ) (67) 

and v(qi,q 3 \ki,k 3 ) is given by 

The amplitude -D3 with the redistribution of colour corresponding to Fig. [5] will be given by 



fdE 1 dE[ fj\ d% d 2 q 2( * , ^ 



Pe-eM^PeA^) < Qu-Qu-Q^qilM^lq'^-q'^-q'^q^ > P E -E[(<£) p E' :L (q'd- ( 68 ) 



The prefactor includes two (— i) from V. 

The amplitude without colour redistribution will differ from D 3a in that between the projectile 
and target pomerons now appear four BKP states 1 1234 >, 1 1243 >, 1 1342 > and 1 1432 > and their 
coupling to the projectile and target pomerons will be the same. This means that now the projectile 
and target pomerons will be connected by 

Mf = \N 2 (v 13 + v 2i - v 23 - «u) [GT i] + Gg 432) + + 4 1342) ] (-13 + -24 - v 23 - v 1A ) 

= \ (H P , 13 + H P , 2A - H P , 23 - H PM ) [G^ + + G^ + Gf^] 



(h PA3 + H P , 24 - H Pj23 - H PM ) . (69) 

The remaining formulas do not change and we find that the amplitude D 3 b without colour redistribu- 
tion corresponding to Fig. [5] will be given by 

ldE 1 dE[ S\j 2/ ' , „ 

Pe-eM^PeA^) < 11, -11,-14^4^^']^, -q[,-q' A ,q' 4 > P E -E[ (q'i)Pe[ il'd- (70) 
To pass to rapidity, consider first evolution of the momentum integrand in time. 



Zl = I dEd ^ 2 e- iEt P E - El P El P E „ E2 P E2 M(E), (71) 



where we suppress the obvious momentum dependence. Presenting 

Pe = j dtP(t)e lEt , M(E) = J dTM(r)e iET 

we have 

Zi(t) = J dhdt 2 dt 3 dUdT J ^^^^-E^+iE^+iiE-E^ts+iE^+iEr 

Pi{t-h)P 2 (t-t 2 )P 3 (t 3 )P 4 (U)M(T)= f dt x /* 1 ^ 2 Pi(t-ti)P 2 (t-ti)M(t 1 -t 2 )P 3 (t2)P4(t2), (72) 



10 Jo 

where we take into account that both P(t) and M(t) are zero at t < 0. Analytically continuing to 
rapidities we find 

ry ryi 

Zi(y) = - / dy x / Pi(y- yi )P 2 (y - y 2 )M( yi - y 2 )P 3 {y 2 )P 4 (y 2 ). (73) 

JO JO 



-1 r\ 



So in the end 



d 3 . = /; *, f *, / ri ^ ( g <; - g « 

-Pi2(y-yi,9i)-P34(y-yi,94) < gi,-?i,-94,?4|MW(yi - y 2 )ki, -9 4 , -9 2 ,9 4 > ^13(2/2, 91)^24(2/2,94) 

(74) 

and 

D » - f *> f * / n ,0) M2 (g ^ - g *: 

-Pi 2 (y - 2/1,91)^34(2/ - 2/1,94) < 91, ~9i, -94,94|M (b) (yi - j/ 2 )|?i, -94, -92, 94 > ^12(2/2, 91)^34(2/2, 94), 

(75) 

where we indicated the numbers of reggeized gluons of which different pomerons are made. 

6 The deuteron-deuteron scattering 

The total deuteron-deuteron scattering cross-section apart from the contributions studied in the pre- 
vious sections will include the standard single and double scattering terms (see Appendix 1. Figs. [TO] 
and fTU) . So the total cross-section is the sum 



a dd = sin 9 le _|_ (jdouble _|_ a (i) _ ^g-j 

i=l 

Here the single cross-section is well known 

a sm9le = (j pp + a nn + 2fj pn . (77) 
The double cross-section is (see Appendix 1.) 

Rouble = _ l _ {(Jpp(Jnn + I d 2 6T 2 (6)) (7g) 

where the transverse density T(b) is expressed in the standard manner via the deuteron wave function: 

T d (b)= [ dz\i; d (b,z)\ 2 . (79) 



Finally the additional cross-sections due to the QCD effects are expressed via Di, i = 1, 2, 3 according 
to Eq. gD 

a« = -2AiV c 2 (<^> d ) 2 , (80) 
where < ... >d means averaging in the deuteron. 

7 Discussion 

For high-energy nucleus-nucleus scattering we have calculated the leading terms in the eikonal function 
for the forward scattering amplitude corresponding to the collision of two scattering centers in the 
projectile nucleus with two scattering centers in the target nucleus. Apart from the obvious pomeron 
exchange, these terms include contributions from connected diagrams involving two pomerons from 
the projectile and two pomerons from the target with all possible interactions in between. We have 
demonstrated that the result is infrared safe, as one expected. It is remarkable that it involves 
contributions from intermediate BKP states formed by 4 reggeized gluons. 

For the deuteron-deuteron scattering the total cross-section is the sum of the eikonal function plus 
the double scattering term in Eq. (|76p . If one takes the pomerons as described by the BFKL equation 



2/ 



Figure 9: Electromagnetic form factor of the deuteron 



then the behavour of the dd cross-section will roughly correspond to the double pomeron exchange 
that is ex.p(2ABFKLy), where Abfkl = (a<jiV c /7r)41n2 is the BFKL intercept. Note that in this case 
the intermediate BKP state in will enter at comparatively low energies, so that one cannot use its 
asymptotical behaviour to find the result. On the other hand description of the total NN cross-section 
by the single pomeron exchange is obviously unrealistic. The single and double cross-sections in @ 
can be calculated using the experimental values of this cross-section. As to the rest, instead of simple 
BFKL pomerons one may use unitarized expressions corresponding to sums of fan diagrams and found 
as solutions of the BK equation. The latter do not grow at high energies and so the asymptotical 
behavour of the cross-section will be determined by that of the BKP state entering a^ 3 \ which grows, 
although not so fast as the pomeron (~ exp 0.2A3AbfklU, [14j). As a result the dd cross-sections 
allow for the direct experimental study of the behaviour of such states, theoretical calculation of which 
presents serious difficulties. 

This is of course true also for collisions of heavy nuclei. But in this case unitarization of the eikonal 
function will in any case lead to the total cross-section which are more or less geometrical. Additional 
terms in the eikonal may of course change it considerably but the cross-section will not change at 
least inside the nucleus where the eikonal remains large. In this case a more interesting problem 
is the inclusive gluon production in the nucleus-nucleus collisions to which disconnected diagrams 
do not contribute. This requires cutting our forward scattering amplitude to select the observed 
intermediate states with due attention to possible cancellation between real and virtual processes 
(AGK cancellations). This problem will be dealt in future studies. 
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9 Appendix 1. Deuteron in the Glauber approach 

9.1 Scattering on the deutron 

To formulate the Glauber approach to the collisions with deuteron in the diagrammatic technique we 
first have to relate the relativistic dpn vertex T with the deuteron wave function. To this end we study 
the electromagnetic form-factor of the deuteron, illustrated in Fig. [91 where vertices T are shown with 
blobs. In the lab. sysstem and at zero momentum transfer it is equal to 2M where M = 2m — e is the 
deuteron mass. So we get the normalization condition 

d% 2Zi r 2 (l 2 ) 

— ,2 Zfl2 ,2 .-nX2 = 2M - ( 81 ) 



((2tt)H (m 2 -q- i0)(m 2 i0)< 




Figure 10: Scattering on the deuteron in the impulse approximation 

Here h + h = 2Z and 4Z 2 = M 2 . We neglect spins and consider all particles as scalar for simplicity. 
We have lo = m — e/2, and 1 = 0, so that putting l 2 = I + A we find 

m 2 — l\ = —2m\o — A^, m 2 — l\ = 2m\o — X±, 

where we used the orders of magnitude Ao ~ e, |Aj_| ~ ^Jme. Integration over Ao transforms (|81 1) into 

d% r 2 (i 2 ) 



(2tt) 3 8m(me + 1 2 ) 2 



2. (82) 



Comparing with the standard normalization of the deuteron wave function ^(1) we find the desired 
relation 

= ,, 9 ,3/2^ 7m . (83) 

which allows to relate the relativistic dpn vertex with the deuteron wave function in the momentum 
space. 

In the impulse approximation, Fig. [lOl with the spectator neutron the corresponding amplitude is 
given by 

Airnp = f J^2 T 2 ^) , . 

J (2vr) 4 i(m 2 -l 2 -i0)(m 2 - Z 2 - i0) 2 ' 1 ' 



where a is the forward scattering amplitude on the proton. Using (|8ip we find that the integral is 
equal to 2 and we get A imp = 2a. But the relativistic flux on the deuteron is twice that on the proton, 
so that we get ad = cr p + a n , where the second term takes into account the diagram of Fig. [TUl with 
the spectator proton. 

Now consider double scattering on the deuteron illustrated in Fig. [TTJ The amplitude is given by 

f d% d% r(i 2 ) rjiy 

J (27r)M27r) 4 i 122 2z \m 2 -q-i0)(m 2 -q-i0)( m i-l[ 2 -i0)( m *-l> 2 2 -i0y 

where H is the high-energy part and it is taken into account that it can only depend on the z- 
component of the transferred momentum, since it is the only of the spatial components which enters 
multiplied by the high projectile momentum. 

Integrations over the zero components of the nuclear momenta factorize and we get 



, d%d% T(h) T(l'2, (m 

A -l-^yr H{hz l2z Um(me + l 2 )4m(me + Y^' m 



or using (|83|) 

A = - I ^fr^H(l 2z - l' 2z )Mh)m'2)- (87) 



1 r d 3 l 2 d% 
mj (2vr) 3 



1 r 



Figure 11: Double scattering on the deuteron 

Integrations over the transverse components are done immediately to convert the wave functions into 
those with the transverse coordinates zero: 

A = - f ^f^H(l 2z - l' 2z )Mr± = 0, hz)Mr± = 0, 4)- (88) 
Transforming completely to the coordinate space we find 

A = - I ^^dzdz'H(q z )Mr± = 0,z)Mr± = 0, z')e u ^ z - z '^ z ' , (89) 

where we introduced the z-component of the transferred momentum putting l' z = l z + q z . Integration 
over l 2z gives our final expression 

A = - I ^dzH(q z )\Mr± = 0,z)\ 2 e~^ z . (90) 
m J 2ir 

The Glauber approximation follows if H(q z ) has a singularity at q z = 0. Typically 

lmH(q z ) = -D(2n)5(2kq) = -D^5(q z ). (91) 

Here we use qo « \q z \ and ko = k z > 0. In this case we get the Glauber approximation for the double 
scattering amplitude 

Im.4 = ~D J dz\Mr± = 0, z)\ 2 = ~D < ^ > d . (92) 
where < ... > d means the average in the deuteron. The cross-section is 

"—h b< 4? >i - (93) 

To see how this formula works consider the simplest case of the double scattering corresponding 
to double elastic collision shown in Fig. [T2l In this case 

H(q z ) = ™P"n( m2 _ {k ~\ z)2 _ i0 + m2 _ {k Z\ z) 2_ J = ^2^(2^), (94) 

so that D = —a p a n . Using ([92]) we find 

ImA = a p a n - < -^—^ > d = -sa p a n < -^—^ > d . (95) 



Figure 12: Double elastic collision on the deuteron 

2k 




Figure 13: Double scattering in deuteron-deuteron collisions 

From this dividing by 2s and doubling to take into account the crossed diagram find the double 
cross-section 

Rouble = _ < 1 >d (96) 

2irr z 

For the nuclear target instead of (|93|) we have at fixed impact parameter b 

a A (b) = -^DTl(b) (97) 
and for double elastic collisions instead of (1961) 



a A (b) = ~A(A-l)a 2 N Tl(b). (98) 



To conclude we note that H and F are both Lorenz invariant. So (|9ip can be used to find F in 
any system. 

9.2 Double scattering in d-d collisions 

Now consider the case when two deuterons collide at high energies and each one experiences double 
collision, illusttated in Fig. [I3j Our treatment is to consider subsequently the two systems in which 
the deuteron is well understandable, the rest systems of the target and projectile deuterons. 
We start from the rest system of the target deuteron. We use Eq. ([92]) and write 



I mA = -- s D 1 <^> d , (99) 



where it is assumed that Im.H\ = —Di2ir5(2kq) and we include into H\ all the rest part of the diagram 
in Fig. [13] including the coupling to the projectile nucleons. 




Figure 14: Double elastic collision in deuteron-deuteron scattering 



Now we boost the system into the rest one for the projectile. In this system we can repeat 
our treatment of the coupling to the two nucleons. If D\ = D2it5(21k) where n is the momentum 
transferred from the projectile, integration over the nucleon momenta will give the same factor (1/s) < 
1/2-Kr 2 >d and we shall get 

W=-(i) 2 £<-L>2. (100) 

Note that we have 

Imil = -D(2ir) 2 5(2(lK))5(2(kq)). (101) 

This means that (jlOOh is symmetric in projectile and target, as expected. From (llOOh we immediately 
get the cross-section (jH) taking in to account the definition of D, Eq. ([2]) 

Special attention is to be given for the case when the high-energy part H is disconnected, shown 
in Fig. [TH Then H contains a 5 function corresponding to conservation laws for the two connected 
parts and is given by 

H = i(a pp a nn + a 2 pn )(2Tr) 4 5 i (k 1 + h - k[ - (102) 
Factor i combines (— i) from the definition of A and i 2 from the two NN amplitudes. Note that 

(27r) 4 <5 4 (fci +h- k[ - l[) = (2k) 4 5(k + + q + )5(K + q„)5 2 (k 1± + l 1± - k' 1± l - l' 1± ) 
= 2s(2ir) 2 5(2(l K ))5(2(kq)) [ d ^b(ki±+h±-k' 1± i-l[ ± ) _ 



As we see, integrations over the transverse coordinates of the nucleons in the projectile and target 
become interdependent. Under the sign of integration over b we include the exponentials depending 
on the transverse momenta of the projectile and target in the corresponding integrals to obtain in 
(JSHj) i)d(b,hz)ipd(b)l'2x) instead of ipd( r ± = 0, l2z)ipd( r ± = (Mzz) an( ^ similarly for the projectile. All 
subsequent calculations remain unchanged and in the end we obtain in (|92p 



i dz\Mb,z)\ 2 = T d (b) (103) 
instead of 

dz\ip d (r ± = 0,z)\ 2 =< ^2 > d ■ 

So the net result of the connection between the transferred transverse momenta is to substitute 
< 1/2-Kr 2 >d—> Td{b) both in the projectile and target and then integrate over b. The rest factors 
from H give D = 2sa p a n and from (llOOj) we conclude 

Im^ = - s {a pp a nn + a 2 pn ) J d 2 bTj(b) (104) 



-i o 



Dividing by 4s we find the cross-section 

a dd = —((TppCTnn + o-pj J d 2 bTj(b) (105) 

which looks very much like the standard Glauber formula. 

For the collision of two heavy nuclei instead of Eq. (|105j) we shall get at fixed b 

<?AB(b) = ~A{A - l)B(B - l)a%( J d 2 b A T A (b A )T B (b - b A )f . (106) 
10 Appendix 2. Colour factors 

The explicit expressions for the colour wave functions of the projectile and target in which pairs 
(12), (34) and (13), (24) respectively form colour singlets are: 

|(12)(34) >= ±5 aia2 5 a3a4 , |(13)(24) >= ±6 aia3 6 a2a4 . (107) 

Here we neglect terms of the relative order 1/N 2 . Their scalar product is 

< (13)(24)|(12)(34) >= -Ls aia3 5 a2a4 5 aia2 5 a;ia4 = ± (108) 

which is the overall damping factor accompanying all diagrams with the redistribution of colour like 
Fig. [21 

We denote C« = —(TjTj). For interactions connecting vacuum pairs either in the projectile or in 
the target CV,- = N c 

< (13)(24)|C ij |(12)(34) >= N c < (13) (24) | (12) (34) >= ^, (ij) = (12), (34), (13), (24) (109) 
For the remaining two interactions we find 
< (13) (24) | Cu | (12) (34) >= ^^ ia ^ ai< <5 aia2< 5 a 3 a4 ^ 02 ^ a 3/<^ c / <a4C = ^1^^ = ~. 

c c c 

(110) 

Interchange (1 f> 2), (3 f> 4) gives 

<(13)(24)|C 23 |(12)(34)>=-i-. (Ill) 

So effectively for these interactions Cy = — N c 

Apart from states | (12) (34) > and | (13) (24) > in the diagrams we encounter six BKP states with 
different ordering of the 4 gluons: 

1 1234 >, 1 1243 >, 1 1324 >, 1342 >, 1 1423 >, 1432 > . 

In the high colour limit their explicit form is 

1 1234 >= -^h aia2C h a3a4C , (112) 

where h abc = d abc + if abc with the properties 

[h ab r = h bac^ Y^h^Yh^ = 5 ab 2N c (l -JL), £ h acd h bcd = -5 ab ^- (113) 

cd c cd c 

The states \ijkl > are cyclic symmetric in (ijkl). 



Their scalar products with the projectile and target states are 



< ( 12) (34) | dj | klmn >= ±~, (116) 



< (12)(34)|1234 >= J^S ai aJa 3 a 4 K ia2C h a3a4C = 0. (114) 

< (13)(24)|1234 >= -^5 aia3 6 a2a4 h aia2C h a3a4C = -L h ^ h a ^ c = -2-L. (115) 

iv c IV c iV c 

Generally if (12) or (34) are neighbors in (ijkl) then states |(12)(fcZ) > and \ijkl > are orthogonal. If 
they are not the the scalar product is the same as in (|115() . 

We also need matrix elements of colour matrices Cjj between projectile (target) states and BKP 
states. Obviously we need only Cij which do not connect vacuum pairs in the projectile (target), 
namely for (ij) = (13), (14), (23), (24) Then we find that for (klmn) = ((1234), (2134), (2143) and 
(1423) 

1 

; 2 ; 

where the sign plus is to be taken when (ij) are neighbours in (klmn) and the sign minus when they 
are not. Acting on the rest two states 1 1324 > and 1 1423 >. all matrices Ci3,C24,Ci4 and C23 give 
N c /2 since neighbour gluons are in the gluon colour state and the matrix elements become damped 
by 1/Nl E.g. 

< (12)(34)|C 13 |1423 >=^< (12)(34)|1423 >= ~. (117) 

(Note that the correct derivation of (|116|) and I117|) in some cases requires taking into account sub- 
dominant terms in (|112p ) 

As a result the matrix elements of C13 are 

< (12)(34)|C 13 |1234 >=< (12)(34)|C 13 |2143 >= ~, 
< (12)(34)|Ci 3 12134 >=< (12) (34) I C 13 1 1243 >= i 



which implies 



< (12) (34) I C13 1 1324 >=< (12) (34) I C13 1 1423 >= 0, 



< (12)(34)|Ci 3 =~\< 1234| - i < 2143| + i < 2134| + i < 1243|. (118) 



Note that the summed probabilities correctly give unity. 
Similarly 

< (12)(34)|C 24 = -- < 1234| - - < 2143| + - < 2134| + - < 1243|, 

< (12)(34)|C 23 = -- < 2134| - - < 1243| + - < 1234| + - < 2143|, 



and 



< (12)(34)|Ci4 = -- < 2134| - - < 1243| + - < 1234| + - < 2143|. 



Interchanging (2 <-> 3) we get 



Ci 2 |(13)(24) >= — -|1324 > — -|3142 > +^|3124 > +i|1342 >, 

C 34 |(13)(24) >= — -|1324 > --|3142 > +^|3124 > +^|1342 >, 
C 23 |(13)(24) >= -i|3124 > -||1342 > +||1324 > +||3142 >, 



Table 1: Matrix elements of the product Cij (lines) by Cki (columns) between states < (12) (34) | and 
| (13) (24) > 





(12) 


(34) 


(23) 


(14) 


(13) 


1/2 





-1/2 


-1/2 


(24) 





-1/2 








(23) 


-1/2 





1/2 


1/2 


(14) 


-1/2 





1/2 


1/2 



C M |(13)(24) >= — - 1 1324 > — -|3142 > +^|3124 > +^|1342 > . 

These relations allow to study matrix elements of the product of two matrices CijCki between the 
projectile and target states. They are shown in Table 1. with lines (ij) and columns (kl) 

From these results we can find the probability to find a particular BKP state between the projectile 
and target. States 1 1234 > and 1 1324 > do not appear and we find the contribution from the double 
interaction VijVki where Vij = —Cijg 2 Vij 

- < (12)(34)|(wi3+U24-V23-t;i4)(|1243 >< 1243+|1342 >< 1342|) (v 12 +v 34 -v 2 3-vu) |(13)(24) > . 

(H9) 

We are also interested in the matrix elements of products of two matrices CijCki between projectile 
and target states without redistribution of color, that is between states < (12) (34) | and | (12) (34) >. In 
particular we shall be interested in separate contribution from BKP states. In this case four different 
BKP states appear between the projectile and target 1 1234 >, 1 1432 >, 1 1342 > and 1 1243 > with 
equal probability and, similar to (|119p we find the probabilities 

- < (12) (34) | ( Ul3 + V 2i - V 23 - V U ) 

(J1234 >< 1234| + |1243 >< 1243+|1432 >< 1432| + |1342 >< 1342|) (v 13 +v 24 -v 23 -vu) |(13)(24) > . 

(120) 
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